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ABSTRACT 

The internal dynamics of a dark matter structure may have the remarkable property that the local 
temperature in the structure depends on direction. This is parametrized by the velocity anisotropy (3 
which must be zero for relaxed collisional structures, but has been shown to be non-zero in numerical 
simulations of dark matter structures. Here we present a method to infer the radial profile of the 
velocity anisotropy of the dark matter halo in a galaxy cluster from X-ray observables of the intra- 
cluster gas. This non-parametric method is based on a universal relation between the dark matter 
temperature and the gas temperature which is confirmed through numerical simulations. We apply 
this method to observational data and we find that is significantly different from zero at intermediate 
radii. Thus we find a strong indication that dark matter is effectively collisionless on the dynamical 
time-scale of clusters, which implies an upper limit on the self-interaction cross-section per unit mass 
<j/m < lcm 2 g _1 . Our results may provide an independent way to determine the stellar mass density 
in the central regions of a relaxed cluster, as well as a test of whether a cluster is in fact relaxed. 

Subject headings: dark matter — galaxies: clusters: general — X-rays: galaxies: clusters 



1. INTRODUCTION 

Our understanding of dark matter structures has in- 
creased significantly over the recent years. This progress 
has mainly been driven by numerical simulations which 
have identified a range of universalities of the dark mat- 
ter structures. One of the first general p roperties to be 
suggested was the radial density profile (|Navarro et al. 
19961 jMoore et al.|[l998| I Diemand et al.||2004| |Merritt 
et al.|2006| Graham et al.|2006| ), whose radial behaviour 
was shown to change from a fairly shallow decline in the 
central region to a much steeper decline in the outer re- 
gions. This behaviour has been confirmed observation- 
ally for galaxy clusters, both through X-ray observations 
dBuote fc Lewis|2004| |Pointecouteau et al.|2005||Arnaud | 
et al.||2005[|Viknlinin et al.||2006[ |Pratt et al.||2006| , and 



zero i n the central region to roughly 0.5 in the outer re 



also more recently through stro ng and weak lensing ob- 
servations dSand et al.|2 004 ; Broadhu rst et al. 2005 Com- 
|erford et~a l. 2006; Limousin et al. 2008p. 

A slightly less intuitive quantity to be considered is the 
dark matter velocity anisotropy defined by 







(1) 



where erf and tr 2 are the 1-dimcnsional tangential and 
radial velocity dispers ions in a spherical system (Bin- 
ney & Tremaine 19871. This anisotropy was shown m 
pure dark matter simulations to increase radially from 
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gion (Carlberg et al. 



1997 Cole & Lacey 1996 Hansen 



& Moore 2006). For collisional systems, in contrast, the 



velocity anisotropy is explicitly zero in the equilibrated 
regions. Therefore, eventually inferring (3 from observa- 
tional data is an important test of whether dark matter 
is in fact collisionless, as assumed in the standard model 
of structure formation. On this note, it has been shown 
that the Galactic velocity anisotropy can affect the detec- 
tion r ates of direct dark matter searches ( Vergados et al 



20081, and it is i n principle measurable in a direction- 
sensitive detector (Host & Hansen 20071. 



The most massive bound structures m the Universe are 
clusters of galaxies, which consist of an extended dark 
matter halo, an X-ray emitting plasma making up the 
intracluster medium (ICM), and the individual galaxies. 
While the contribution of galaxies to the total mass is 
small, approximately 10 % of the cluster mass resides 
in the ICM. The present generation of X-ray satellites, 
XMM-Newton and Chandra, allows very accurate mea- 
surements of azimuthally-averaged radial profiles of den- 
sity and temperature of the ICM. These are used, under 
the assumption of hydrostatic equilibrium and spherical 
symmetry of both gas and total mass distributio ns, to 



estimate total, ga s, and dark matter mass profiles (Fab 



ricant et al]|1980 1 

Below we infer the radial velocity anisotropy profile of 
dark matter in 16 galaxy clusters using a generally ap- 
plicable framework without any parametrized modeling 
of the clusters. In short, we assume a universal relation 
between the effective temperature of dark matter and the 
ICM temperature, which allows us to solve the dynamics 
of the dark matter halo using the radial gas temperature 
and density profiles determined from X-ray data. We 
investigate the shape and validity of this temperature 
relation in two cosmological simulations of galaxy clus- 
ters, based on independent numerical codes. We apply 
our method to 16 galaxy clusters from two different sam- 
ples and find a velocity anisotropy significantly different 



2 



from zero in the outer parts, in qualitative agreement 
with simulations. 

Our approach here is a generalization of th e non- 
parametric analysis in Hansen & Piffaretti ( 20071 where 
(3 was inferred neglecting the radi al dependence. We also 
note the parametr ized analyses in |Ikebe et al.| ( |2004[ ) and 



0, or 



Cole & 



Morandi & Ettori ( 2007 ), where the total dark matter ve 
locity dispersion was infer red assuming eithe r /3 
the analytica l /3-profil es of|Colm et al.|(|2000[ ) or 
|Lacey| (119961) (see also|Wojtak et al.| ( |2008| )). In particu 
lar, |Morandi fc Ettorij p007[ ) found that the dark matter 
temperature and the ICM temperature were essentially 
the same in strong cooling-core clusters. 

The structure of the paper is the following: In the next 
section, we discuss how we relate the temperature of dark 
matter to the observable gas temperature. In section 3 
we show how we can then solve the dynamics of the dark 
matter. In section 4 we test the assumed temperature 
relation and our method on numerical simulations, and 
in section 5 we apply the method to observational data. 
Section 6 is the summary and discussion. 

2. THE TEMPERATURE OF DARK MATTER 

The equality of inertial and gravitational mass implies 
that the orbit of a test particle in a gravitational sys- 
tem is independent of mass. For example, the veloc- 
ity of a circular orbit in a spherical mass distribution 
v 2 c = GM(r) jr depends only on the distance to the center 
of the system and the mass contained within that radius. 
Therefore it is natural to assume that, at a given radius, 
all species in a relaxed, spherical gravitational system 
have the same average specific kinetic energy. Obviously, 
they also have the same specific potential energy. In a gas 
system, equilibrium implies energy equipartition between 
all species. It is clear that the corresponding principle 
for a relaxed gravitational system is a common velocity 
dispersion, precisely because gravitational dynamics are 
independent of mass. Since the average velocity is asso- 
ciated with the thermal energy content, this relationship 
is expressed by 

Tdm = KT gas . (2) 

The parameter k is constant as long as the impact of 
radiative or entropy-changing processes affecting the gas 
is negligible and the system is relaxed. Therefore, we al- 
low for a radial dependence, k — K(r/r25oo), where r25oo 
is the scale radius within which the mean total density 
is 2500 times the critical density at the redshift of the 
cluster. 

The dark matter temperature in ^ is naturally not 
well-defined as there is no thermodynamic equilibrium 
for a collisionless gas. Instead, we simply define an ef- 
fective dark matter temperature which is proportional to 
the three-dimensional velocity dispersion, 



= o/ im » (°"r + 2cr t 



(3) 
(4) 



The velocity dispersion has been decomposed into the 
contributions from the one-dimensional radial and tan- 
gential dispersions. We choose the constant of propor- 
tionality to be the mean molecular mass of the intraclus- 
ter gas simply to allow n to be of order unity. Equations 



(pjl-Q are equivalent to assuming that the specific ener- 
gies of gas and dark matter particles are the same up to 
a factor of n, on average. The same conjecture was made 



Hansen & Piffaretti ( 2007 1 but with n = 1 explicitly, 
t should be mentioned that the temper ature relation 



( 2 ) wa s recently analyzed in simulations by Evrard et al. 
(7008). Whereas we allow a possible radial variation in 
the temperature relation, those authors considered aver- 
ages within r 2 oo and found that 



fcsT gas //xmij 



b <''20Q 



= 1.04 ±0.06, 



(5) 



'DM 



This was based on their determination of n 
(0.87 ± 0.04) (r spcc /T mw ), 
troscopic temperature to 



<'' , 2()0 



where the ratio of the spec- 
the mas s-weighted t empera- 
1.1 ± 0.05 ( |Nagai|[2006l ). Then, 
Rincs et al. (|2008|) it was noted that by appl ying 



tu re was (r sP ec/ T; 



virial scaling to t he WMAP5+SN+BAO results (Ko 



matsu et al.| 2008[ ), an average value of K _1 | <r . B00 = 1.1 
was found. The authors concluded that the observational 
results indicated that the average specific energy of the 
ICM was close to both that of the dark matter and that 
of the galaxies. In section [4] we will arrive at the same 
conclusion for simulated galaxy clusters. 

3. SOLVING THE DARK MATTER DYNAMICS 

Equation (p| allows us to estimate the total velocity 
dispersion profile of the dark matter structure from mea- 
surements of the radial temperature profile of the gas. In 
this section we discuss how we can proceed to determine 
the dark matter velocity anisotropy. 

The collisionless Jeans equations relate the dynamical 
properties of the dark matter to the gravitational poten- 
tial of the cluster. Assuming that the system is spheri- 
cally symmetric and in a steady st ate, the second Jeans 
equat ion can be put in the form (Binney & Tremaine 
1987| 



dr 



2vl 



GM 



(6) 



where v is the dark matter number density, vf is the 
second moment of the ith velocity component, and M 
is the mass contained within radius r. If it is further 
assumed that there are no bulk flows, ul = 0, and that the 
tangential velocity dispersions are equal, a\ = a^ = of, 
the Jeans equation becomes 



dlnp DM 
din r 



din of 
din r 



2/3)=- 



GM(r) 



(7) 



where pdm is the mass density, a% is the radial velocity 
dispersion, and /? is the velocity anisotropy introduced 

in 0-. 

Similar to the Jeans equation, the radial part of the 
Euler equations of the ICM expresses the condition that 
the thermal pressure of the gas balances the gravitational 
potential. This equation of hydrostatic equilibrium reads 



din n e 
din r 



dlnTi 



ga.s 



dlnr 



GM(r) 



(8) 



where T gas is the gas electron temperature and n e is the 
number density of electrons. This important equation 
has been widely used to estimate the total mass of a 



3 



galaxy cluster from X-ray data. In case there is turbu- 
lence or larger scale bulk motion in the gas a dditional 
terms of the form (v ■ V)v — (vg + v^) jr appear (Landau 
fc Lifshitz||1987 1. Neglecting such terms may lead to an 
underestimate of the mass, however this is usu ally not a 
major effect for sy stems that appear relaxed ( |Piffaretti 
fc Valdarn ini 2008). 

By equating ^fy and p| and using ([I]) and ([2| to elim- 
inate (3, we obtain the following differential equation for 
the radial velocity dispersion, 



0V 



rflnflpM 
d In r 



d In of 
d In r 



+ 3 =ip(r), 



(9) 



where the function tp is defined by 

ip(r) = 3k . 

[imn r 

Clearly, ip is determined directly from the X-ray observ- 



(10) 



ables and the K-profile, which we discuss in section 4.3 

The differential equation (9]) is solved by finding an 
integrating factor which yields 



pT>u(r) r 3 J 



dr' ip{r')p BU (r'y 



12 



(11) 



The dark matter density is determined as usual through 
Pdm = Ptot — nmHfi e . With the radial velocity disper- 
sion profile determined, the velocity anisotropy is easily 
recovered from either the temperature relation Q or the 
Jeans equation §7^, 



(12) 
(13) 



2/5 tr = 31 



2/3jo 



dln(pBM<rr) 



din r 



GM 



Obviously these two expressions should be equal. This 
can be used as a consistency check on whether numeri- 
cal issues related to the differentiations and integration 
involved are kept under control. 

To summarize, the assumed relation (p| between the 
effective dark matter temperature and trie gas temper- 
ature, along with the mass estimate from (TSJ> , allows us 
to solve the dark matter dynamics directly from X-ray 
data, and determine both the radial velocity dispersion 
and the velocity anisotropy as functions of radius. 

4. CLUSTER SIMULATIONS 

We use numerical simulations of the formation of 
galaxy clusters in the ACDM cosmology to investigate 
the validity and shape of the temperature relation ([2]), 
and to test the method for determining the velocity an- 
isotropy In order to check systematic effects we take 
samples from two different simulations based on two com- 
pletely independent numerical codes. 

4.1. CLEF 

We first consider a sample of 67 clu sters taken from 
the CLEF simulation (|Kay et al.||2007[) , details of which 
are briefly summarized here. The CLEF simulation was 



run w ith the GADGET2 N-body/SPH code (Springel 



|2005J ) and followed the evolution of large-scale structure 
within a box of comoving length, 200/i _1 Mpc. The fol- 
lowing cosmological parameters were assumed: f2 m = 



0. 3; n A = 0.7; fi b = 0.0486; h = 0.7; n s = l;cr 8 = 0.9. 
Here the value of the Hubble constant is written as 
100 /ikms" 1 Mpc -1 and erg is the rms mass fluctuation 
at the present epoch in a sphere of radius 8/i _1 Mpc. 
The number of particles was set to 428 3 for each of the 
gas and dark matter species, thus determining the par- 
ticle masses to be ttidm = 7.1 x 10 9 /i _1 Mq and m gas = 
1.4 x 1Q 9 /i _1 M q respectively. The equivalent Plummer 
softening length was set to 20/i _1 kpc and held fixed 
at all times in comoving co-ordinates. Pressure forces 
were calculated using the standard GADGET2 entropy- 
conserving version of SPH with an artificial viscosity to 
convert kinetic energy into thermal energy where the flow 
was convergent. The gas could cool radiatively assuming 
a fixed metallicity, Z = 0.3Z Q . Cold (T < 10 5 K) gas 
with Tin > 10~ 3 cm -3 either formed stars or was heated 
by an entropy, AS = lOOOkeVcm 2 . This choice was de- 
termined stochastically by selecting a random number, r, 
from the unit interval and heating the particle if r < 0.1, 

1. e. a 10 per cent probability of being heated. This high 
level of feedback was necessary to reproduce the o bserved 
excess entropy in clusters (see Kay et al. (20071 for fur- 
ther details). 

To select the cluster sample, we first consider all clus- 
ters at z = with X-ray temperatures, kT > 2keV; 
this produces 95 objects, with virial masses, M vir > 
1.3 x 10 14 /i _1 M Q (correspondingly, > 15,000 dark mat- 
ter particles) . We then select those clusters with 3D sub- 
structure statistic, s < 0.05. The substructure statistic 



(Thomas et al. |1998| measures the displacement of the 
centre of mass from the potential minimum of the cluster 
(taken to be its centre), relative to r^oo, which is the scale 
radius within which the mean total density is 500 times 
the critical density. By making this cut, we therefore ex- 
clude all clusters that show significant signs of dynamical 
activity, i.e. major mergers. 

4.2. V06 

The second s ample is a subsample of the one presented 



Valdarnini (2006) which we refer to as V06. These 



simulations assumed a concordance flat ACDM with the 
same cosmological parameters as for the CLEF simula- 
tion. 

The simulation ensemble of galaxy clusters was con- 
structed according to a procedure described in ( Piff aretti 
fc Valdarnini|2008| . Here we briefly summarize the most 
important aspects. The hydrodynamic simulations were 
run using an entropy-conserving multistep TREESPH 
code for a sample of 153 clusters spanning a range from 
~ 1.5 x 10 15 h- 1 M G down to M mr ~ 1.5 x lO 14 ^ 1 ^. 
The initial conditions (z;„ = 49) were extracted from a 
set of purely N-body cosmological simulations in which 
clusters of galaxies were identified from the particle dis- 
tribution at z = using a friends-of-friends algorithm. 
In order to investigate the effect of the implemented gas 
processes on the energy equipartition between gas end 
dark matter particles, we performed both adiabatic and 
radiative simulations. The radiative simulations are of 
course more realistic than the adiabatic ones, because 
they additionally take into account radi ative cooling, sta r 
formation, energy and metal feedback ( |Valdarnini|2003[ | . 
More details concerning the simulation technique and the 
i mplementation of p hysical processes of the gas are given 
in |Valdarnini| ( |2006[ ). 
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Fig. 1. — Radial profile of K = Tbivi/Tgas for the samples of 
clusters obtained from the CLEF and V06 simulations compris- 
ing 67 and 20 clusters, respectively. We plot the median and la 
percentiles taken over each sample. The vertical line indicates the 
largest radius of the observational data sample, while the vertical 
lines indicate the mean and standard deviation of the K-profile that 
we use in the fiducial analysis. Note that, for the CLEF sample, 
only eight clusters contribute to the innermost bin. 

In order to avoid contamination from dynamically per- 
turbed clusters, we select the 20 most relaxed objects at 
z = 0. The selection is based on the power ratio method, 
which measures the amount of substructure in X-ray 
surface brightness maps. The map sources a pseudo- 
potential which is expanded in plane harmonics, and the 
ratio of the third coefficient to the zeroth i s a measure 
of substructure. More details are given in Piffaretti fc 



Valdarnini (20081 



4.3. 



The temperature relation 

We examine the temperature relation ^ in the two 
simulated samples by comparing the gas mass-weighted 
temperature to the rescaled dark matter velocity disper- 
sion. The resulting K-profiles are shown in fig. [I] and 
clearly n sa 1 for both samples. Since we apply some- 
what different criteria to select the two simulation sam- 
ples, it is not surprising to find slightly different pro- 
files. This indicates a systematic uncertainty of ±0.1 in 
the simulated k profiles. The kinetic energy associated 
with bulk motions of both gas and dark matter particles 
is at most a few percent of the thermal energy within 
2 r25oo i outside which bulk motion is n ot negligible. Thi s 
is in agreement with what was found in |Ascasibar (2003 1. 
Due to the standard problem of limited force resolution, 
the simulations do not probe the innermost region reli- 
ably. Therefore we exclude data inside a cutoff radius 
(56/i _1 kpc for CLEF, 0.1r 2 5oo for V06), which means 
we cannot estimate k in the central region where gas 
physics can make a significant impact. 

The adiabatic version of the V06 sample exhibits 
a larger median K-profile which is constant about 1.2 
within r25oo and increases steadily t o 1.4 at r2no- This i s 
comparable with the earlier work of Rasia et al. (2004 1, 
where the specific energy of dark matter was seen to be 
larger than that of the gas by 20-30% in adiabatic sim- 
ulations. 



4.4. Reconstructing the velocity anisotropy 

In order to test the method outlined above for deter- 
mining (3, we reconstruct the anisotropy profiles observed 
in the simulated samples. Here, we assume k = 1 for all 
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Fig. 2. — Comparison of estimated and true values of the physical 
quantities involved in determining the velocity anisotropy f3 in our 
simulations. Top, the ratio of the reconstructed total density to 
the true; bottom, the ratio of the reconstructed a^, to the true one. 
Error bars show the la percentiles taken over the sample members. 



radii even though we expect deviations at small radii. 
First we derive the integrated mass profile M(r) for each 
cluster assuming hydrostatic equilibrium d8J), and from 
that the total density profile. The numerical derivatives 
involved are calculated using three-point quadratic inter- 
polation. The estimated density profile, shown in fig. [2] 
(top), displays a satisfactory agreement with the actual 
density profile for both the CLEF and V06 samples. The 
only exception is at the outermost radii above r 2 5oo where 
the density is underestima ted . Next, we calculate the 
radial velocity dispersion (111 by interpolating the in- 
tegrand from r = using a four-point natural spline 
interpolation. We compare the resulting radial velocity 



dispersion with the actual in fig. [2] (bottom) which shows 
that there is good agreement except for the deviation at 
large radii already seen in the density profiles. 

Finally we determine the velocity anisotropy parame- 
ter (3. We find similar results whether we calculate /3j or 
Ptn however the temperature relation yields less noisy re- 
sults. The median velocity anisotropy profiles are shown 
in fig. [3] together with the median actual profile. The 
reconstructed profile tracks the actual anisotropy well in 
the inner parts but overestimates (3 in the outer parts. 
There is also considerable noise in the results. 

In order to understand the origin of the deviations at 
large radii and the significant scatter in our results, we 
investigate the systematics of the analysis, as applied to 
the CLEF sample (similar conditions hold for the V06 
sample). First, we substitute the dark matter density es- 
timated from hydrostatic equilibrium with the true den- 
sity. The /3-profiles calculated on this basis are shown in 
the top panels of fig. |4j The agreement between the esti- 
mated and actual (3 is considerably improved, and the er- 
ror bars are significantly reduced. This clearly indicates 
that, in the fiducial analysis, the numerical derivatives 
necessary to estimate /3dm are responsible for the large 
error bars. Since we do not want to do any parametrized 
modeling of the gas properties, the numerical derivatives 
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Fig. 3. — Reconstructed velocity anisotropics for the simulated 
samples. The hatched bands show the actual /3-profiles of the sam- 
ples. Error bars show the ltr percentiles taken over the sample 
members. 



are liable to amplify noise and induce systematic devia- 
tions in the outermost bin, where the quantities are only 
constrained to one side. Additionally, this explains why 
/3j e appears more noisy in the fiducial analysis since an 
additional derivative must be calculated. The test also 
shows that there is a deviation from hydrostatic equi- 
librium at large radii which is part of the reason why (3 
is overestimated. As a second test, we additionally use 
the true three-dimensional velocity dispersion instead of 
using the temperature relation. This yields further im- 
provement as to how well the reconstructed /3 tracks the 
true one, as shown in the bottom panels of fig. [4j This 
implies that it is possible to get the correct scale of the 
radial velocity dispersion, calculated as an integral from 
the center, despite the lack of resolved data in the inner 
radii. We note that, with respect to observational data, 
the tests we apply here can possibly be utilized in the 
future, e.g. with accurate density profiles inferred from 
gravitational lcnsing, and with more detailed knowledge 
of k from improved simulations. We conclude that the 
numerical simulations provide proof that our method is 
robust and that it is indeed possible to infer the /3-profile 
despite lacking knowledge of k in the center. 

5. OBSERVATIONS 

Next we apply our analysis to observational data from 
which the radial gas density and temperature profiles 
are recovered. This is done strictly using non-parametric 
methods, i.e. no modeling of the gas properties is in- 
volved. Our data consists of the deprojected density and 
temperature profiles of two samples of clusters at low and 
intermediate redshift, respectively. The deprojected pro- 
files were obtained from X-ray data analysis published in 
earlier work (details below). We consider clusters which 
appear relaxed and close to spherical, and for which suf- 
ficient spectroscopic data are available to analyze several 
annuli, so that the radial variations of the gas density 
and temperature are resolved with good statistics. 

The first set of eleven clusters at low redshift is 
based on X-ray data from XMM-Newton of the clusters: 
A262, A496, A1795, A1837, A2052, A4059, Sersic 159-3, 
MKW3s, MKW9, NGC533, and 2A0335+096. These ob- 



Fig. 4. — Systematics of the reconstruction of the /3 p rofi les for 
the CLE F simulation. Again, (3 is recovered both from ( | 1 2[ i (left) 
and l |13| (right). The true dark matter density is substituted for the 
estimated, and in the bottom panels we additionally use the true 
total velocity dispersion instead of estimating it from Tdm = KT gas . 



jects are highly relaxed cool-core (CC) clusters selected 
as to match the requirements described above. The ob- 
j ects w ere part of the sample analyzed in Kaastra et al. 
(2004} (see this paper for an extensive presentation of the 
data analysis), in which deprojected radial temperature 
and density profiles were derived from spatially resolved 
spectroscopy . We adopt the radial bin selection of |Pif-| 



faretti et al. (20051 in order to ensure a robust determi- 
nation of gas temperature and density for the full radial 
range. Note that data for A2052 and Sersic 159—3 were 



also used in the analysis by Hansen & Piffaretti (20071 



where a constant velocity amsotropy was assumed 

The other set of five intermediate redshift X-ray galaxy 
clusters (RXJ1347.5, A1689, A2218, A1914, A611) is 



from the Chandra sample analyzed in Morandi et al 
(2007). The radial deprojected temperature and density 



profiles were retrieved through resolved spectral analy- 
sis in a set of annuli, selected to collect at least 2000 net 
counts, by assuming spherical geo metry and by using the 
definition of 'effective volume' (see Morandi et al. (20071 
for further details). 

6. RESULTS 

We determine the dark matter velocity anisotropy pro- 
file j3{r) of each cluster according to the recipe in section 
[3] using a Monte Carlo method. For each radial bin the 
deprojected gas temperature and density are sampled as- 
suming Gaussian uncertainties, i.e. a random number is 
chosen from a Gaussian distribution with mean equal to 
the estimated temperature or density and a standard de- 
viation equal to the uncertainty of the estimate. The 
bins are sampled independently. The parameter k is also 
sampled for each bin, assuming a Gaussian distribution 
with a mean of 1 and a standard deviation of 0.1, which 
is a reasonable value according to the simulations. The 
sampled profiles are used to reconstruct the total mass 
through (8j, and then the integrand, the radial velocity 
dispersion, and the velocity anisotropy are calculated in 
each bin. The sampled set of profiles is accepted only if 
the temperature and density as well as the reconstructed 
dark matter density and radial velocity dispersion are all 
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Fig. 5. — Three steps in the calculation of the velocity anisotropy for Sersic 159—3. Left, the inferred total density; center, the radial 
velocity dispersion; right, the /3-profiles. The gas density and temperature profiles are also shown. The scale radius for this cluster is 
estimated to be r25oo = 337 ± 13 kpc. Error bars indicate the propagated statistical uncertainties on the ICM temperature and density 
profile, taken as the la percentiles of 1000 Monte Carlo samples. This is unlike in the previous figures where the error bars indicate the 
spread over the numerically simulated samples. In the right panel, the radial positions of /3tr and /3j e have been offset slightly for clarity. 



non-negative in all bins. For each sample, we also esti- 
mate the scale radius ^2500 and the mass M2500 contained 
within that radius. Table [T] summarizes the properties of 
the clusters in our sample. 

The numerical methods for calculating derivatives and 
integrals are the same as for the simulated samples, 
i.e. three-point quadratic interpolation is used for deriva- 
tives and f our- point spline interpolation is used for the 
integral in ( flT) . The integration results are stable to us- 
ing two-point linear, three-point quadratic, or four-point 
least squares quadratic interpolation instead. 

Individual steps of the reconstruction are shown in 
fig.[5]for the cluster Sersic 159—3, and the deprojected in- 
put data are also displayed. We always plot the median 
and la percentiles since spurious outliers in individual 
Monte Carlo samples can bias the mean and standard de- 
viation significantly. The size of the error bars is mostly 
determined by the uncertainties of the temperatures, to 
a lesser degree by the uncertainties of the ICM densi- 
ties, and it is virtually insensitive to the 10% variation 
assumed for the K-profile. 

As can be seen in the right panel of fig. [5] the agree- 
ment between /3 t r and /3j G indicates that numerical effects 
associated with the integration and differentiations are 
small. On the other hand, /3 becomes unphysically large 
in the outermost bins since the reconstructed radial ve- 
locity dispersion for some samples becomes greater than 
the total velocity dispersion. This result is similar to 
that found in the blind analysis of the simulation sam- 
ples. As discussed above, this behaviour is mainly due 
to a deviation from hydrostatic equilibrium of the gas, 
and to a lesser degree to edge effects making the numeri- 
cal differentiations less well determined in the outermost 
bin. It is possible that systematic uncertainties in the in- 
put data or radial variations in k for individual clusters 
also play a role. In principle, we could impose a 2 < a 2 , 
thereby forcing j3 < 1, as another physical condition on 
each Monte Carlo sample, but we prefer not to do so in 
order to have a consistency check. 

We repeat the data analysis for the remaining 15 clus- 
ters of our sample and the resulting velocity anisotropy 
profiles are shown in fig. [6] In almost all cases the an- 
isotropy is small in the inner radial bins and increases to 
between 0.5 and 1.0 in the outer parts. There is good 
agreement between the two derivations of (3 for all clus- 
ters, indicating that numerical issues are under control. 



Since the qualitative behaviour of the velocity aniso- 
tropy profiles are similar, we combine all our data into 
a single 'stacked' profile, shown in fig. [7] In the re- 
gion where direct comparison is possible, the measured 
stacked profile is very similar to the reconstructed f3 
profiles for the simulation samples (the green line), and 
within r25oo there is also agreement with the actual veloc- 
ity anisotropy of the simulation samples (hatched band) . 
The velocity anisotropy is likely overestimated outside 
?*2500 f° r the same reason as for the simulated samples, 
i.e. deviation from hydrostatic equilibrium, but the ef- 
fect appears to be even stronger for the observational 
data. Interior to the cut-off radius of the numerical 
simulations, the observations tend to (3 ~ 0.3. This is 
somewhat surprising since numerical simulations at all 
mass scales generally have very little anisotropy towards 
the center of structures. While we cannot exclude the 
possibility that cluster halos are anisotropic even at low 
radii, our result can also be explained by the neglected 
stellar contribution p± to the total mass density. To 
first order, this contribution enters our analysis in the 
Jeans equation through the estimated dark matter den- 
sity pdm = Pdm + P*- In terms of 5* = p±/puM, the 
Jeans equation becomes 



din Pdm 
dlnr 



d In a 2 
din r 



2/3- 



c2m(l + <y 
dlnr 

GM{r) 



(14) 



where the slope of (1 + 6+) is negative since the stellar 
density must fall off faster than the dark matter density. 
This means that we overestimate the velocity anisotropy 
in the central region by not accounting for the stellar 
mass. Indeed, if we assume that 50% of the total mass in 
the innermost bin is made up of stars, the velocity aniso- 
tropy in the two innermost bins becomes consistent with 
zero. There is also a seco nd order correction through 
the appearance of pdm in (11) instead of pdm , but this 



correction must be small since the density contributes to 
both the integrand and the normalization factor. 

Finally, we investigate how the assumed shape of the k- 
profilc affects our results. We try five different profiles as 
functions of x = r/r25oo with noise added as before, and 
calculate the velocity anisotropy profiles for each. The 
K-profilcs are chosen so as to mimick either the effects 
of gas radiative cooling or AGN heating in the central 



7 



1.5 

1.0 

0.5 

0.0 

-0.5 
1.5 

1.0 
0.5 
0.0 

-0.5 
1.5 

1.0 
0.5 
0.0 



RXJ1 347.5 

r 25 oo= 707 kpc 



ll 



A1689 

r z5 oo= 612 kpc 



i HI 



1 




A1914 



591 kpc 



tlttfl 



A611 



520 kpc 



A1837 



374 kpc 



t \ 




A262 

<2soo= 255 kpc 



♦ill' 



A4059 

hsoo = 445 kpc 



MKW9 

r 2500 = 280 kpc 



.* til 



1 



A1795 

r 2500 = 504 kpc 



\ 



A2052 

1-2500 = 362 kpc 

Mil' 


: MKW3S 

|-25oo = 403 kpc 


: NGC533 

|-25oo= 191 kpc 

I 


ill- 


i i ; 


1 1 I ' I 1 ■ 1 


i i i ' i i i ■ i 



1.00 0.01 



0.10 

^ r 2500 



1.00 0.01 



0.10 



1.00 



Fig. 6. — Median velocity anisotropy profiles for the remaining 15 clusters of our sample. The estimated scale radii are also shown, and 
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Fig. 7. — Median velocity anisotropy profile of all 16 clusters in 
our dataset. In this case the error bars denote the la percentiles 
of the combined probability density of all clusters within the bin. 
The actual and reconstructed /3-profiles from the simulations are 
also shown. The left vertical line is the innermost radius probed in 
the CLEF simulations and the right vertical line shows, roughly, 
the onset of significant deviations from hydrostatic equilibrium in 
the simulations, see fig. [2] 



TABLE 1 
Properties of our cluster sample 



Cluster 



^2500 /kpc 



M 2 tmo/M q 



A262 


0.015 


256 ± 


28 


(2.7 ± 


0.8) x 


10 13 


A496 


0.032 


398 ± 


10 


(1.0 ± 


0.2) x 


10 14 


A1795 


0.064 


504 ± 


22 


(1.9 ± 


0.2) x 


10 14 


A1837 


0.071 


374 ± 


26 


(8.0 ± 


1.7) x 


10 13 


A2052 


0.036 


362 ± 


11 


(6.7 ± 


0.6) x 


10 13 


A4059 


0.047 


445 ± 


21 


(1.3 ± 


0.2) x 


10 14 


Sersic 159-3 


0.057 


337 ± 


17 


(5.7 ± 


0.8) x 


10 13 


MKW3s 


0.046 


404 ± 


14 


(9.5 ± 


0.9) x 


10 13 


MKW9 


0.040 


279 ± 


44 


(3.2 ± 


1.5) x 


10 13 


NGC533 


0.018 


191 ± 


15 


(9.7 ± 


2.2) x 


10 12 


2A0335+096 


0.034 


350 ± 


40 


(6.9 ± 


2.5) x 


10 13 


A611 


0.29 


519 ± 


52 


(2.5 ± 


0.6) x 


10 14 


A1689 


0.18 


609 ± 


4 


(3.5 ± 


0.7) x 


10 14 


A1914 


0.17 


590 ± 


44 


(3.3 ± 


0.8) x 


10 14 


A2218 


0.18 


535 ± 


51 


(2.5 ± 


0.7) x 


10 14 


RXJ1347.5-1145 


0.45 


710 ± 


60 


(7.3 ± 


1.4) x 


10 14 



regions, or to check the results if the dark matter is gen- 
erally hotter or cooler than the gas. The radially varying 
profiles we try are extreme cases of the simulation pro- 
files, fig. [T] Typically, the result is that the /3-profile 
is shifted in the central regions while the outer regions 
are largely unaffected, as shown in fig. [8} This analysis 
confirms that there is a significant velocity anisotropy 
at large radii, independent of the specific assumptions 
about the temperature relation. 

7. SUMMARY AND DISCUSSION 

In this paper, we have presented a non-parametric 
method to infer the velocity anisotropy of dark matter in 
clusters of galxies from the observable temperature and 
density of the intracluster medium. We assume that the 
intracluster medium has the same specific energy as the 
dark matter, and we investigate the validity of this as- 
sumption in two different cosmological simulations of the 
formation of galaxy clusters. Both confirm the simplest 
possible form of the relation, namely 7dm ~ T sas in the 
radial range which is resolved. 

We have tested how well our method can reconstruct 



the actual velocity anisotropy in the simulated clusters, 
and we have found good agreement between the two, 
although the reconstruction is sensitive to systematic bi- 
ases connected with deviations from hydrostatic equilib- 
rium. 

We have applied our method to the radial ICM density 
and temperature profiles of 16 galaxy clusters based on 
Chandra and XMM-Newton X-ray data. The shape of 
the velocity anisotropy profiles is always consistent with 
that seen in simulations, which tends to zero at the inner- 
most radius where the temperature relation is calibrated. 
It then increases to about 0.5 at r25oo and even larger in 
the outer regions. The same is true of the fiducial anal- 
ysis applied to simulated data and is likely caused by a 
deviation from hydrostatic equilibrium outside T25oo- We 
also find a significant anisotropy even if we assume radi- 
ally varying K-profiles, such as can be expected given the 
strong gas cooling and AGN heating in the core of many 
clusters, or if we assume k ^ 1. The agreement between 
the observed velocity anisotropy and that predicted in 
numerical simulations shows that we are beginning to 
understand also the dynamical aspects of dark matter in 
halos. 

In the innermost radial bins we measure a rather large 
anisotropy, but this is most likely an overestimation due 
to the neglect of the stellar mass in the center. This can 
be used as a means to estimate the stellar mass profile 
of galaxy clusters if one assumes that the velocity dis- 
persion to be isotropic in the central regions. Similarly, 
our method may be used as a general test of whether a 
cluster is relaxed. A reconstructed velocity anisotropy 
which deviates significantly from the simulated profiles 
would be a strong hint that the data do not support the 
assumption of hydrostatic equilibrium. 

The inferred velocity anisotropy profiles are signifi- 
cantly different from zero which means that the collec- 
tive behaviour of dark matter is unlike that of baryonic 
particles in gases. This shows that dark matter is ef- 
fectively collisionlcss on the timescale of r ~ 10 9 , the 
dynamical timescale of galaxy clusters. By taking typ- 
ical values at ~ 0.3r 2 5oo and allowing only a few scat- 
terings within the time r, this corresponds to an order- 
of-magnitude upper limit to the scattering cross-section 
of roughly ajm — ((OrjjVlTu) -1 ^ lcm 2 g _1 . This limit 
i s similar to what has been found for mer ging clusters 
( |Markevitch et al |2004| iBradac et aL|2008| ), and within 
an order of magnitude of the scattering cross-section for 
self-in t eracti ng dark matter proposed in |Spergel fc Steim] 



hardt| ( |2000 l. 

We emphasize that improvements to the numerical 
simulations in the near future will improve our under- 
standing of the k profile and hopefully track the impact 
of radiative effects in the center. We also hope that 
improved understanding of deviations from hydrostatic 
equilibrium will allow us to estimate how large the sus- 
pected bias at large radii is. On the observational side, 
the main problem at present is the uncertainty in the 
temperature profile. Improvements can be expected both 
with regards to the deprojection analysis and the amount 
of data available. Obviously, there is also the possibility 
of including a kincmatical analysis of the galaxy clusters 
in our method. 
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Fig. 8. — The effect of assuming different K-profiles on the stacked velocity anisotropy profile. Top left, the five K-profiles. Others, the 
resulting sample averaged /3-profiles calculated assuming the numbered K-profile. In all cases, j3 is greater than zero in the outer parts. 
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